How to resum long-distance contributions to the QCD pressure? 
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The strict coupling constant expansion for the free energy of hot QCD plasma shows bad conver- 
gence at all reasonable temperatures, and does not agree well with its 4d lattice determination. This 
has recently lead to various refined resummations, whereby the agreement with the lattice result 
should improve, at the cost of a loss of a formal agreement with the coupling constant expansion 
and particularly with its large infrared sensitive "long-distance" contributions. We show here how 
to resum the dominant long-distance effects by using a 3d effective field theory, and determine their 
magnitude by simple lattice Monte Carlo simulations. 
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Introduction. At temperatures above 200 MeV, the 
properties of matter described by the laws of QCD are 
expected to change. The system should look more like 
a collection of free quarks and gluons than a collection 
of their bound states, such as mesons. It is a chal- 
lenge to find observables which would clearly manifest 
this change, and hopefully also be directly or indirectly 
measurable in heavy ion collision experiments. 

From the theoretical point of view, one of the simplest 
observables witnessing the change is the free energy of 
the plasma, or its pressure 0. Indeed, according to the 
Stefan-Boltzmann law, the value of the free energy counts 
the number of light elementary excitations in the plasma, 
be they quarks and gluons, or mesons. 

The reality is somewhat more complicated. Interac- 
tions change the Stefan-Boltzmann law, so that pressure 
is no longer proportional to the number of degrees of free- 
dom. And in fact, interactions are strong. An explicit 
computation of the free energy to order 0{g^T'^) j^-^] 
shows that there are large corrections, with alternating 
signs, such that convergence is poor at any reasonable 
temperature. Of course, at least without light dynami- 
cal fermions, the full pressure can still be obtained with 
4d finite temperature lattice simulations However, 
in order to really understand the properties of the QCD 
plasma phase, one would also like to have some analytic 
understanding of the origin of this result. 

A way of at least understanding why the convergence 
is poor is the observation that when as = g'^/{4Tr) ^ 
1, the system undergoes dimensional reduction 
and its static long wavelength "soft" or "light" degrees 
of freedom can be described by a three-dimensional (3d) 
effective field theory, 

Csd = ^Tr + Tr [A, ^0]' + m|,Tr + AA(Tr A^)^ 
where ~ g^T^,XA ~ g'^T are parameters com- 



puted perturbatively up to optimised next-to-leading or- 
der level (see below). This effective theory is confining, 
therefore non-perturbative ]lO| , pT| . In Q C^d was used 
to reproduce the perturbative free energy up to order 
0{g^T^) and the bad convergence was shown to be 
due precisely to these degrees of freedom. 

Our objective here is to study the free energy of QCD 
by including the dominant badly convergent contribu- 
tions from Csd non-perturbatively, to all orders, by using 
lattice Monte Carlo simulations. In this way, we can find 
out how important the combined effect of the badly con- 
vergent series really is in the free energy. 

It is important to keep in mind that infrared sensi- 
tive effects can be different in various quantities. For 
instance, the free energy is dominated by ultraviolet de- 
grees of freedom, and the long-distance effects we study 
here may turn out to be subdominant. Thus it would 
be wrong to conclude that any approach which manages 
to reproduce the numerical data for the free energy in 
a satisfactory way, would also reproduce other quanti- 
ties. A good testing ground for this are the longest static 
correlation lengths in the QCD plasma: they are fully 
non-perturbative, but it is already known that the re- 
sults of 4d simulations jl2j are reproduced precisely by 
the infrared degrees of freedom we employ in C^d I^,||,|l3| ■ 

The relation of our approach to the other recent ap- 
proaches for the determination of the free energy of 
QCD p^p^ can be described as follows. At present, 
these approaches do not reproduce the known 0{g^T'^) 
result in the limit of a weak coupling, nor do they account 
for any genuine non-perturbative contributions. Thus 
large infrared effects are suppressed without an a priori 
justification; the justification comes a posteriori through 
the reasonable agreement with numerical data. Our re- 
sults here attempt at providing a theoretical understand- 
ing for why the long-distance contributions need not be 
important in the QCD pressure. 



1 



Method. The pressure or the free energy density 
of QCD is a quantity which formaUy gets contributions 
both from short-distance physics {I 1^{t^T)^^) and long- 
distance physics (/ ^{gT)~^). The separation of the free 
energy into these two different types of contributions has 
been discussed in detail in Q. Interactions between the 
short and long-distance modes account for the param- 
eters of the effective long-distance theory C^^, and in 
addition there is an additive part coming directly from 
the short-distance modes, as we will presently specify. 

To describe the effects of the short-distance modes in 
detail, we find it useful to introduce the dimensionless pa- 
rameters y = mj-j /g^,x = XA/g^, where g'^ is the gauge 
coupling within the effective theory. In terms of the phys- 
ical parameters T,Ay^ of QCD, next-to-leading order 
"fastest apparent convergence" optimised perturbation 
theory tells that |^ (for a number of flavours Nf = and 
colours Nc = 3), 
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The result of Eq. (36), can now be expressed as 
follows. Using the MS scheme with the scale parameter 
/i3d, let us compute the dimensionless quantity 
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where V is the volume. The pressure can then be ex- 
pressed as (we have here again put Nf = 0, iV^ = 3) 
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where po{T) = {n^T'^ /AbjiN^ - 1 + {7/4)NcNf) is 
the non- interacting Stefan-Boltzmann result. The /isd- 
dependence here is cancelled by that in !F^{x,y). 

A few comments on Eq. (^ are in order. First of 
all, the term proportional to y could also be written as 
^ 0{x'^), and at the present level of accuracy there is 
no unique way of making a distinction. We have cho- 
sen the present form because the relatively large loga- 
rithmic term is then dealt with in connection with ^j^, 
whereby cancellations occur. Second, strictly speaking 
ln(^3d/r) should be replaced with ln(/23d/T) -|- 6, but 
S = 7£;-ln2-41/2160-(17/72) ln27r-(37/36)[lnC]'(2)-f 
(19/72)[lnC]'(4) w 1.35x10-4 can be ignored for all prac- 
tical purposes. Finally, with the expressions available 
at present, the relation in Eq. (^ has an error start- 
ing at order 0{g^), corresponding to 0(1/(477)^) within 
the parentheses. This correction is however from short- 
distance physics alone, and we shall ignore it here. 
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FIG. 1. The pressure in Eq. (^, with the long-distance part 
from Eq. 1^ included to various loop orders. The 4d lattice 
results are from the first reference in Q. It should be noted 
that they have a normalisation ambiguity at low temperatures 
T <Tc allowing for a small shift of the curve. 



By Eqs. @), (||), (|j), the perturbative short-distance 
contribution to the pressure has been accounted for to 
a satisfactory level, and we are left with evaluating the 
long-distance part, T-^{x,y). The perturbative expres- 
sion for T-^{x, y) is known up to 3-loop level, corre- 
sponding to 0(5^T4) accuracy in piT). Adding terms 
involving the scalar self-interaction x to the result of [Q , 
we can write 
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where dA — — 1,Ca — Nc, and AT-^{x,y) accounts 
for the higher order corrections. In terms of the 4d cou- 
pling constant, all contributions involving x in Eq. ^ are 
of order 0{g^) or higher, while the terms y^ ,y\ny, y^ 
are of orders 5'^, Infl/g), g^, respectively. 

As is well known the convergence of the pertur- 

bative expansion in Eq. (|^) is quite poor when values of 
X, y corresponding to any reasonable physical tempera- 
ture r/A^|g are chosen. For future reference, we illus- 
trate this in Fig. |^. We have used Eqs. 
together with terms up to order y 2 in Eq. (^ 

The idea of our approach of improving the determina- 
tion of Ty^{x,y) is the following. We write 



ATy^{x,y) = AJ'j;jg(xo,?;o) 
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where y — y{x) is defined in Eq. (Q). The partial deriva- 
tives are now given by adjoint Higgs field condensates: 

"y \ 53 / MS \ 93 I MS, pert 

where A^j g^y^ ^^^^ is the perturbative result up to 

0(y~^), obtained by taking a derivative of Eq. (||) with 
respect to y. In the case of d ^Ty^j dx , a similar relation 
is obtained but with the condensate ((Tr Aq)^). 

On the other hand, with a computation in lattice 
perturbation theory, a condensate measured in lattice 
Monte Carlo simulations can be related to the con- 
densates (TrylQ)j;jg, ((TrA§)^)^. Due to the super- 
renormalizable nature of /^sd , such analytic relations can 
be computed exactly near the continuum limit |p7| , |l8| . 

Thus, we need to evaluate the condensates on the lat- 
tice, transform the result to the MS scheme, and per- 
form finally the integration in Eq. (|^) numerically. When 
added to A^^(a;o, 2/0)7 we obtain a non-perturbative re- 
sult, which we can plug into Eq. (^. 

What remains is to determine the integration constant 
A^j^(a;o, j/o)- The idea is that despite the bad conver- 
gence shown in Fig. |l], at high enough temperatures the 
form of AjFj^(a;o, ?/o) is known. Indeed, inspecting the 
general structure of Eq. (||) , we know that 

A^^(:ro,yo) = 7f^dACl(l + 0(^,-^)). (8) 

Here eo, containing an unknown logarithmic dependence 
on 2/0, represents the famous non- perturbative 0{g^T^) 
term Suppose now that we choose T = To ~ 

lO^^A^, corresponding to = 1.0 x 10~^, yo = 3.86. 
Then the higher order terms in Eq. are expected to 
be subdominant, since Cyi/(47ryo^) ~ 0.1 and xq/Ca ^ 
0.01, and we only need to know eo- 

The main error sources of this non-perturbative and 
unambiguous setup are as follows: 

(a) Even though in principle an independent non- 
perturbative determination of eo is possible for instance 
by measuring the condensate (Tr Ff^ ) along the lines 
in |l9[ , doing this systematically requires a 4-loop compu- 
tation in lattice perturbation theory, and this is beyond 
our scope here. Therefore we will treat eo as a free inte- 
gration constant whose magnitude will be fixed below. 

(b) Due to the smallness of x/Ca: we will also ignore 
here the term arising from dAT^^/dx in Eq. (^). 

(c) The numerical procedure introduces small statisti- 
cal errors, as well as systematic errors from the extrapo- 
lations to the infinite volume and continuum limits. 

(d) Finally, we should of course remember that the 
effective theory C^d loses its accuracy when higher order 
operators not included become important. In fact, for 
Nf = the QCD phase transition is related to the so 
called Z(3)-symmetry [^^, and this symmetry is not 
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FIG. 2. The diflterence in Eq. (g). Here /3g = d/igh), 
where a is the lattice spacing, and the continuum limit corre- 
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FIG. 3. The pressure after the inclusion of A!F-^{x, y) 
from Eq. (^. Statistical errors are shown only for eo = 10. 

fully reproduced by [p|j2l|], without all the higher 
order operators. There are many indications, however, 
that the effective theory should be rather accurate down 
to low temperatures, T ~ 2Tc [^|J|,|l^. Below that, some 
other effective description may apply (see, e.g., 

Numerical results. After this background, we show 
in Fig. H the difference in Eq. (0), measured with lattice 
simulations. This result is then used in Eq. (H) to obtain 
ATy^{x,y). When added to Eqs. (|), (|), we obtain 
Fig. As discussed above, the boundary value at (al- 
most) infinite temperature, determined by eo, is for the 
moment a free parameter. 

We observe that at low temperatures the outcome de- 
pends strongly on the value of eo- The correct value 
would appear to be eo ~ 10.0 ± 2.0. Even then, the 
present results lose their accuracy at T ~ 5Tc, but seem 
to work well above this. Exploiting the full power of the 
dimensionally reduced theory down to its limit T ~ 2Tc 
would also necessitate the inclusion of ((TrAg)^). 
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Discussion. In 4d lattice simulations, there is a (nu- 
merically small) ambiguity in the determination of the 
pressure, because only pressure differences can be mea- 
sured, and thus an integration constant has to be spec- 
ified at low temperatures in a non-perturbative regime. 
Here we fix the integration constant by starting from the 
opposite direction, from very high temperatures. This 
allows to determine all quantities in terms of T/ and 
the number of fermion flavours, without ambiguities. We 
can also address a huge range of temperatures, unlike 4d 
simulations which can only go up to T a few x Tc- 

The result of our procedure is summarised by Eqs. (^)- 
and Fig. |^. We draw two important conclusions. The 
first is that the outcome depends strongly on the non- 
perturbative contribution of order 0(g^T'^) Jo[, as can 
be observed from the eo-dependence in Fig. |. The value 
of Co could in principle be determined by a well-defined 
procedure, although in practice it is a project of con- 
siderable technical complication. But our present study 
provides an estimate for what the result should be. The 
order of magnitude 0(10) seems reasonable, since it is 
known from other contexts such as the Debye mass fl^ ] 
that non-perturbative constants tend to be large. 

The second is that when the large non-perturbative 
0{g^T'^) term is summed together with the set of all 
higher order terms determined via (TtAq), then these 
long-distance contributions almost cancel at T > 30Aj^! 
Indeed, the sum, the curve with cq ~ 10 in Fig. | , does 
not differ there much from the term 0{y^) in Fig. l[ For 
smaller temperatures, 1^:,T < 30Aj^, on the other 

hand, only our numerical results are trustworthy. 

Finally, we also find that although the dependence on 
the effective scalar self-coupling x is of high perturbative 
order, in practice it is expected to play a role as one 
approaches Tc. Its contribution can be obtained from the 
condensate ((Tr Ag)^). To relate this to the MS scheme 
requires again a perturbative 4-loop computation. 

Let us end with a philosophical note. When one wants 
to understand 4d simulation results, one could argue 
that one should aim at almost fully analytical resum- 
mations ||l^-[l^. However, we suspect that these are un- 
avoidably specific for the particular observable consid- 
ered: they may work for the entropy or pressure because 
the result is short-distance dominated, but would fail for 
instance for Debye screening where long-distance effects 
are dominant. It seems to us that it may ultimately 
be more useful to obtain a unified understanding of the 
relevant degrees of freedom in the system, even if some 
observables have to be evaluated numerically. 
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